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Abstract 

A posteriori estimates for mixed finite element discretizations of the 
Navier-Stokes equations are derived. We show that the task of estimating 
the error in the evolutionary Navier-Stokes equations can be reduced to 
the estimation of the error in a steady Stokes problem. As a consequence, 
any available procedure to estimate the error in a Stokes problem can be 
used to estimate the error in the nonlinear evolutionary problem. A prac- 
tical procedure to estimate the error based on the so-called postprocessed 
approximation is also considered. Both the semidiscrete (in space) and 
the fully discrete cases are analyzed. Some numerical experiments are 
provided. 

1 Introduction 

We consider the incompressible Navier-Stokes equations 

u t -Au+(u- V)it + Vp = /, (1) 
div(u) = 0, 

in a bounded domain SI C M. d (d — 2, 3) with a smooth boundary subject to 
homogeneous Dirichlet boundary conditions u = on 9fi. In (pj), u is the 
velocity field, p the pressure, and / a given force field. For simplicity in the 
exposition we assume, as in [5], [27], [25], [2J5], [33J, that the fluid density and 
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viscosity have been normalized by an adequate change of scale in space and 
time. 

Let Uh and ph be the semi-discrete (in space) mixed finite element (MFE) 
approximations to the velocity u and pressure p, respectively, solution of (TTJ) 
corresponding to a given initial condition 

u(.,0) = uo. (2) 

We study the a posteriori error estimation of these approximations in the L 2 
and H 1 norm for the velocity and in the L 2 /R norm for the pressure. To do this 
for a given time t* > 0, we consider the solution (u, p) of the Stokes problem 

-Au + Vp = f- j- t u h {t*) - (u h (t*) ■ VH(t*)l ^ ^ 

div(w) = 0j m ' (3) 

u = 0, on d£l. 

We prove that u and p are approximations to u and p whose errors decay by 
a factor of /i|log(/i)| faster than those of Uh and ph (h being the mesh size). 
As a consequence, the quantities u — Uh and p — Ph, ar e asymptotically exact 
indicators of the errors u — Uh and p — ph in the Navier-Stokes problem (JIJ— ((2J). 

Furthermore, the key observation in the present paper is that (uh,Ph) is 
also the MFE approximation to the solution (u,p) of the Stokes problem ©. 
Consequently, any available procedure to a posteriori estimate the errors in a 
Stokes problem can be used to estimate the errors u — and p — Ph which, as 
mentioned above, coincide asymptotically with the errors u—Uh and p~ph in the 
evolutionary NS equations. Many references address the question of estimating 
the error in a Stokes problem, see for example [2] , [6] , [7] , [32] , [35] , [39] , [40] and 
the references therein. In this paper we prove that any efficient or asymptotically 
exact estimator of the error in the MFE approximation {Uh,Ph) to the solution 
of the steady Stokes problem ([3]) is also an efficient or asymptotically exact 
estimator, respectively, of the error in the MFE approximation (uh,Ph) to the 
solution of the evolutionary Navier-Stokes equations ([I])-©. 

For the analysis in the present paper we do not assume to have more than 
second-order spatial derivatives bounded in L 2 (fl) d up to initial time t = 0, 
since demanding further regularity requires the data to satisfy nonlocal com- 
patibility conditions unlikely to be fulfilled in practical situations [27] , [28] . The 
analysis of the errors u — u and p — p follows closely [16] where MFE approxi- 
mations to the Stokes problem ([3|) (the so-called postprocessed approximations) 
are considered with the aim of getting improved approximations to the solution 
of |T])-([2]) at any fixed time t* > 0. In this paper we will also refer to (u, p) 
as postprocessed approximations although they are of course not computable in 
practice and they are only considered for the analysis of a posteriori error esti- 
mators. The postprocessed approximations to the Navier-Stokes equations were 
first developed for spectral methods in [53], [2J], [TB], [SS] and also developed 
for MFE methods for the Navier-Stokes equations in [4], [5], [16] . 

For the sake of completeness, in the present paper we also analyze the use 
of the computable postprocessed approximations of [16j for a posteriori error 
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estimation. The use of this kind of postprocessing technique to get a posteriori 
error estimations has been studied in [19] , [20] and [15] for nonlinear parabolic 
equations excluding the Navier-Stokes equations. We refer also to [33] where 
the so-called Stokes reconstruction is used to a posteriori estimate the errors 
of the semi-discrete in space approximations to a linear time-dependent Stokes 
problem. We remark that the Stokes reconstruction of [33J is exactly the post- 
processing approximation (u,p) in the particular case of a linear model. 

In the second part of the paper we consider a posteriori error estimations 
for the fully discrete MFE approximations U% « Uh{t n ) and P£ ~ ph(t n ), 
(t n = At„_i for n — 1,2,..., N) obtained by integrating in time with either 
the backward Eulcr method or the two-step backward differentiation formula 
(BDF). For this purpose, we define a Stokes problem similar to ([3]) but with 
the right-hand-side depending now on the fully discrete MFE approximation U£ 
(problem (fT0"]) - (jTT]) in Section@]below). We will call time-discrete postprocessed 
approximation to the solution (U n , P n ) of this new Stokes problem. As before, 
(U n , P n ) is not computable in practice and it is only considered for the analysis 
of a posteriori error estimation. 

Observe that in the fully discrete case (which is the case in actual com- 
putations) the task of estimating the the error u(t n ) — UV; of the MFE ap- 
proximation becomes more difficult due to the presence of time discretization 
errors — Uh(t n ) — UJf, which are added to the spatial discretization errors 
u{t n ) — Uh(t n ). However we show in Section 3] that if temporal and spatial er- 
rors are not very different in size, the quantity U n — U£ correctly esimates the 
spatial error because the leading terms of the temporal errors in U n and U% get 
canceled out when subtracting U n — U%, leaving only the spatial component of 
the error. This is a very convenient property that allows to use independent 
procedures for the tasks of estimating the errors of the spatial and temporal dis- 
cretizations. We remark that the temporal error can be routinely controlled by 
resorting to well-known ordinary differential equations techniques. Analogous 
results were obtained in [15] for fully discrete finite element approximations to 
evolutionary convection-reaction-diffusion equations using the backward Euler 
method. 

As in the semidiscrcte key point in our results is again the fact 

that the fully discrete MFE approximation (Ufi, P£) to the Navier-Stokes prob- 
lem (JT])— ([2]) is also the MFE approximation to the solution (U n ,P n ) of the 
Stokes problem (l70l) - (|7T]) . As a consequence, we can use again any available 
error estimator for the Stokes problem to estimate the spatial error of the fully 
discrete MFE approximations (U£,Pfi) to the Navier-Stokes problem ([I])-©. 

Computable mixed finite element approximations to (U n ,P n ), the so-called 
fully discrete postprocessed approximations, were studied and analyzed in [17] 
where we proved that the fully discrete postprocessed approximations maintain 
the increased spatial accuracy of the semi-discrete approximations. The anal- 
ysis in the second part of the present paper borrows in part from [17] . Also, 
we propose a computable error estimator based on the fully discrete postpro- 
cessed approximation of [17] and show that it also has the excellent property of 
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separating spatial and temporal errors. 

The rest of the paper is as follows. In Section 2 we introduce some prelimi- 
naries and notation. In Section 3 we study the a posteriori error estimation of 
semi-discrete in space MFE approximations. In Section 4 we study a posteri- 
ori error estimates for fully discrete approximations. Finally, some numerical 
experiments are shown in Section 5. 



2 Preliminaries and notations 

We will assume that fl is a bounded domain in M. d , d = 2,3, of class C ,n , for 
m > 2. When dealing with linear elements (r = 2 below) ft may also be a 
convex polygonal or polyhedral domain. We consider the Hilbert spaces 

H = {u G L 2 (0) d I div(it) = 0, u ■ n Ul = 0} , 
V = {u G Hl{0) d | div(u) = 0} , 

endowed with the inner product of L 2 (0) d and _ffo(0) d , respectively. For I > 
integer and 1 < q < oo, we consider the standard spaces, W Lq (Q) d , of functions 
with derivatives up to order I in L q {tt), and H l (Vt) d = W la {Vl) d . We will 
denote by || • ||; the norm in H l (£l) d , and || • ||_; will represent the norm of its 
dual space. We consider also the quotient spaces H (f2)/K with norm \\p\\w /r — 
inf{||p + c||j | ce R}. 

We recall the following Sobolev's imbeddings [I]: For q G [1, oo), there exists 
a constant C — C(f2, q) such that 

\\v\\ Lq - <C\\v\\ W s.^ ^>--^>0, q<oo, veW s ' q {Sl) d , (4) 
q q a 

For q 1 = oo, @ holds with i < |. 

The following inf-sup condition is satisfied (see [3S]), there exists a constant 
> such that 

ml sup , > p, (5) 

gei2(n)/m„ ei? i (fi)< i ||v||i||?||x, 3 /]k 

where, here and in the sequel, (•, ■) denotes the standard inner product in L 2 (fl) 
or in L 2 {fl) d . 

Let II : L 2 (0) d — > H be the L 2 (0) d projector onto H. We denote by A 
the Stokes operator on 0: 

A:V{A)<ZH — > H, A = -IIA, 2?(A) = H 2 (0) d n V. 

Applying Leray's projector II to (JT|), the equations can be written in the form 

u t + Au + B(u,u) =Uf in 0, 

where B(u, v) = H(u ■ V)v for u, v in Hq (ft) d . 
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We shall use the trilinear form b(-, •, •) defined by 

b(u, v, w) = (F(u, v), to) Vtt, v, w G (f2) d , 

where 

F(u, w) = (u • V)t> + i(V • «)t> Vm,!)£ Hi (Q) d . 

It is straightforward to verify that b enjoys skew-symmetry: 

b(u, v, w) = -b(u, w, v) yu,v,w e Hl(fl) d . (6) 

Let us observe that B(u,v) = UF(u,v) for u G V, v G i?o(fi) d . 

Let us consider for a G K and t > the operators A Q and e — , which are 
defined by means of the spectral properties of A (see, e.g., [13j p. 33], [2Tj). 
Notice that ^4 is a positive self- adjoint operator with compact resolvent in H. 
An easy calculation shows that 

\\A a e- tA \\ <(ae- 1 ) a t- a , a > 0, t > 0, (7) 

where, here and in what follows, ||-|| when applied to an operator denotes the 
associated operator norm. 

We shall assume that the solution u of CD)-© satisfies 

||u(t)||i<Mi, ||u(t)|| 2 <M 2 , 0<<<T, (8) 

for some constants Mi and M 2 . We shall also assume that there exists a constant 
M 2 such that 

||/||i + ||/t||i + ||/tt||i <M 2 , 0<t<T. (9) 
Finally, we shall assume that for some k > 2 

L(fe-2)/2j 

sup ||5 t Lfe/2J /|| fc _ 1 _ 2 , fc/2 | + V] sup ||3 t J /|| < +0O, 

o<*<t" "* 1 2L ' £/2J rr( o<t<T M "* 2j 2 

— j=o — 

so that, according to Theorems 2.4 and 2.5 in [27], there exist positive constants 
M k and K k such that the following bounds hold: 

||u(t)|U + ||«t(*)||fc-2 + lb(i)IU-7« ^ MferW 1 -^ 2 , (10) 

/ <T k -3(s)(\\u(s)\\ 2 k + \\u s (s)\\l_ 2 + |b(s)||^_ 1/R + \\p s (s)\\ 2 Hk - 3/R ) ds < Kl (11) 
Jo 

where r(t) = min(i, 1) and a n = e~ a ^*~ s - ) r n (s) for some a > 0. Observe that 
for t < T < oo, we can take r(i) = i and er„(s) = s". For simplicity, we will 
take these values of r and a n . 

Let Tk = (t/ 1 , ^ )iei h , h > be a family of partitions of suitable domains 
f2/t, where /i is the maximum diameter of the elements r/ 1 G 7^, and $ are the 
mappings of the reference simplex tq onto t,- . 
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Let r > 2, we consider the finite-element spaces 

Sh,r = {xhtC (& h ) | X h\r t O $ G P r_1 (T0)} C ff 1 ^), Sl r = S^H-ffo (0fc), 

where P r_1 (ro) denotes the space of polynomials of degree at most r — 1 on 
ro. As it is customary in the analysis of finite-element methods for the Navier- 
Stokes equations (see e. g., jS], [27], [25], [25], [3D]) we restrict ourselves to 
quasiuniform and regular meshes Tk, so that as a consequence of |12[ Theorem 
3.2.6], the following inverse inequality holds for each Vh € (Si r ) d 

\\vh\\w™><*(n h y <Ch l ~ m ~ d (~~*)\\v h \\ wl , q , {nh)d , (12) 

where 0<l<m<l,l<q'<q<oo. 

We shall denote by (Xh, r , Qh,r—x) the so-called Hood-Taylor element [51151]. 
when r > 3, where 

Xh,r = (S%, r ) d , Qh,r-i = S h , r . 1 nL 2 {0 h )/R, r > 3, 

and the so-called mini-element [10] when r = 2, where Qfc,i = Sft,2 H L 2 (0/i)/K, 
and Xh,2 — 2 ) d ®B/j. Here, B^ is spanned by the bubble functions b T ,r £Th, 
defined by b T (x) = (d + l) d+1 Ai(a;) ■ • • Xd+i(x), if x £ r and elsewhere, where 
Ai(x), . . . , Xd+i(x) denote the barycentric coordinates of x. For these elements 
a uniform inf-sup condition is satisfied (see [9j), that is, there exists a constant 
/3 > independent of the mesh grid size h such that 

. (q h , V • w h ) 

ml sup — j|— J, — j| > p. (13) 

vheQh,r-i Vh ex Kr II u /i||i||*i||l 2 /ir 

We remark that our analysis can also be applied to other pairs of LBB-stable 
mixed finite elements (see Qj)l Remark 2.1]). 

The approximate velocity belongs to the discrete divergence-free space 

Vh,r = X Kr n { Xh G Hl(@ h ) d I (q h , V • X h) = Vq h G Q V -l} . 

which is not a subspace of V. We shall frequently write Vh instead of Vh, r 
whenever the value of r plays no particular role. 

Let : L 2 (0) d — > Vh, r be the discrete Leray's projection defined by 

(U h u,Xh) = (u,Xh) Vxh€Vh, r . 
We will use the following well-known bounds 

||(/-n h Hli < Ch l -i\\u\\ h i</<2, j = o, l. (14) 

We will denote by Ah '■ Vh — > Vh the discrete Stokes operator defined by 
(Vv h , V<f> h ) = (A h v h , <t> h ) = (Al /2 v h , Al /2 h ) Vv h , <t> h G V h . 
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Let (u,p) G (H 2 (n) d DV) x (if^OyR) be the solution of a Stokes problem with 
right-hand side g, we will denote by s/, = S/i(ii) € Vh the so-called Stokes pro- 
jection (see [35]) defined as the velocity component of solution of the following 
Stokes problem: find (sh,qh) G (Xh,r, Qh,r-i) such that 

(Vs h , V<j) h ) + (Vq h , 4> h ) = (g, 4> h ) V(t> h G X Kr , (15) 
(V ■ S h ,lP h ) = (16) 

The following bound holds for 2 < £ < r: 

||«-*fc||o + ft||«-a/i||i < C/i'(||u||i + ||p|| ff! -i /R ). (17) 

The proof of (fTTj) for £1 = can be found in [28 . For the general case, 
flh must be such that the value of 5(h) = max xe gQ h dist(x,dfl) satisfies 5(h) — 
0(h 2 ^ r ~ 1 ^). This can be achieved if, for example, dfl is piecewise of class C 2 ' 1 * -1 -*, 
and superparametric approximation at the boundary is used [3] . Under the same 
conditions, the bound for the pressure is [25] 

\\p - Qkhym < C^h l - l (\\u\\i + Ibll^.-i/a), (18) 

where the constant Cp depends on the constant j3 in the inf-sup condition (|13[) . 
We will assume that the domain f2 is of class C m , with m > r so that standard 
bounds for the Stokes problem [3], [32] imply that 

WA^UgW^. <\\g\\., -l<j<m-2. (19) 

For a domain fl of class C 2 we also have the bound (see [11] ) 

\\p\\m/R<c\\g\\ . (20) 

In what follows we will apply the above estimates to the particular case in 
which (u,p) is the solution of the Navier-Stokes problem (d])-©. ^ n ^at case 
Sh = Sh(u) is the discrete velocity in problem (|15[) — (pT6|) with g = f— u t — (u-Vtt). 
Note that the temporal variable t appears here merely as a parameter, and then, 
taking the time derivative, the error bound (|17[) can also be applied to the time 
derivative of Sh changing u, p by u t , Vt- 

Since we are assuming that f2 is of class C m and m > 2, from (TlTl) and 
standard bounds for the Stokes problem [3j [22] , we deduce that 

\\(A- 1 U- Ay^Uh) f\\. < Ch 2 ~ J \\f\\ V/eL 2 ((]f, j =0,1. (21) 

We consider the semi-discrete finite-element approximation (uh,Ph) to (u,p), 
solution of (P)-©. That is, given u h (0) = U h u , we compute u h (t) £ X hr and 
Ph(t) G Qh,r-i, t G (0,T], satisfying 

(iih,4>h) + (Vu h ,V4>h) +b(uh,Uh,4>h) + (Vpft, cj>h) = (f,<fih) ^ 4>h G X h , r , 

(22) 

(V • uh, iph) = V4e Qh, r -i- 

(23) 
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For 2 < r < 5, provided that (HTJ) - (HHJ) hold for I < r, and (fTD ]) -([TT ]) hold 
for k = r, then we have 

l|u(*)-«fc(*)Ho + /i|Kt)-Ufc(*)||i<C^ ?F 5572 , 0<i<T, (24) 
(see, e.g., [M CHI US]), and also, 

lbft)-^ft)lUv«^ C ' f (r fe -2)/2 - 0<i<T, (25) 

where r' = r if r < 4 and r 1 = r + 1 if r = 5. 
see [29] Proposition 3.2]. 

3 A posteriori error estimations. Semidiscrete 
case 

Let us consider the MFE approximation (uh,Ph) at any time t* £ (0,T] to 
(u(t*),p(t*)) obtained by solving (|2"2")l - (1231) . We consider the postprocessed ap- 
proximation (u(t*),p(t*)) in (V, L 2 (il)/R) which is the solution of the following 
Stokes problem written in weak form 

(Vu(i*),V0) + (Vp(f),0) = (/, 0) - 6(«h(t*), « fc (t*), ^) - ^I?6) 

(V-«(t*),^) = 0, (27) 

for all £ i?o(^) d and V € L 2 (fl)/R. We remark that the MFE approximation 
(u h (t*),p h (t*)) to (u(t*),p(t*)) is also the MFE approximation to the solution 
(tt(t*),p(t*)) of the Stokes problem fl2H])-([I7]). In Theorems [Q and H below 
we prove that the postprocessed approximation (u(t*),p(t*)) is an improved 
approximation to the solution (u,p) of the evolutionary Navier-Stokes equations 
{T])-((2]) at time t*. Although, as it is obvious, (u(t*),p(t*)) is not computable 
in practice, it is however a useful tool to provide a posteriori error estimates 
for the MFE approximation (uh,Ph) at any desired time t* > 0. In Theorem [1] 
we obtain the error bounds for the velocity and in Theorem [2] the bounds for 
the pressure. The improvement is achieved in the H 1 (0) d norm when using the 
mini-element (r = 2) and in both the L 2 (£l) d and i? 1 (J7) d norms in the cases 
r = 3,4. 

In the sequel we will use that for a forcing term satisfying © there exists 
a constant M 3 > 0, depending only on M 2 , \\AhUh(0)\\o and sup 0<t<T 
such that the following bound hold for < t < T: 

\\A h u h {t)\\ 2 < Ml (28) 

The following inequalities hold for all Vh,Wh € Vh and <f> £ -ffg(f2) d , see [29j 
(3.7)]: 

\b(v hl v hl 4>)\ < C ||^||? /2 ||A^||J /2 ||0||o, (29) 

\b(vh,Wh,4>)\ + \b(wh,Vh,4>)\ < c\\vh\\i\\AhW h \\o\\<t>\\a- (30) 
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The proof of Theorem Q] requires some previous results which we now state 
and prove. 

We will use the fact that ||vlj^ 2 M;/,|| = ||VtOft|| for Wh € Vh, from where it 
follows that 

C-^A-^VHo < IKIU < C||A- 1/2 ^|| Vw h eV h , (31) 
where the constant C is independent of h. 

Lemma 1 Let (u,p) be the solution of (QP~fiP and fix a > 0. Then there exists 
a positive constant C = C(M2,a) such that for w^w"^ G Vh satisfying the 
threshold condition 

\\w l h -u\\j <ah 3 / 2 -\ j =0,1, 1 = 1,2, (32) 

the following inequalities hold for j = 0, 1; 

\\A- j/2 n h (F(wlwl)-F(wl,wl))\\ < C\\At i)/2 «-™l)\\ , (33) 
\\A- j/2 R h (F(wlwl) - F(u,u))\\ Q < Cll^-itH^.. (34) 

Proof Due to the equivalence (f3"Tj) and and since ||II/i/||o < ||/||o for / € 
L 2 (il) d it is sufficient to prove 

WFiwlwD-Fiw^W^KCWwl-wh-j, j = 0,1, (35) 

for w = w\ or w — u. We follow the proof [5J Lemma 3.1] where a different 
threshold assumption is assumed. We do this for w — w\, since the case w = u 
is similar but yet simpler. We write 

F(wlwl) ~ F{wl,w 2 h ) = F(wle h ) + F(e h ,w 2 h ), (36) 

where e^ = w\ — w\ . We first observe that 

H^e^wDllo = sup 

H0llo = l 

< C\\e h \\ L 2 d \\Wwl\\ L 2 d /( d -i) +C||e^||i||u;, 2 l ||Lo= 

L 2d/(d-i) + ||ui 2 ||l=) ||e?t||i) 

where, in the last inequality, we have used that thanks to Sobolev's inequality 
(jH) we have ||e/ l || I; 2<j < C7]|e^|| x- Similarly, 

\\ F ( w L e h)\\o < C||*4||z,°°||eft||i +C||Vw / 1 l || i 2d/(d-i)||e/ l || L 2 £i 

< C(\\w h \\ L ~ + WVwlW^nt-^WehWi. 

The proof of the case j = in f|35|) is finished if we show that for I = 1,2, 
both ||w/J|l°° and || Vu>^|| £ 2d/(d-i) are bounded in terms of M-i and the value a in 



(e h -Vu£),0) + -((V-e h )u#),0) 
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the threshold assumption (j3"2")l . To do this, we will use the inverse inequality (|12p 
and the fact that the Stokes projection Sh = Sh(u) satisfies that 

IMU~<C a , \\Vs h \\ L 2« <C S 

for some constant C s = C S (M2) (see for example the proof of Lemma 3.1 in [5]). 
We have 

lKlk~ < \\w l h - a h \\ L oa + \\s h \\ L ~ < Ch- d ' 2 \\w l h - flfcHo + || S/l || L », 
where in the last inequality we have applied (|12l) . and, similarly, 

||V^|| L 2d/(d-l) < \\V{w l h - S h )\\ L 2d/(d-l) +\\Vs h \\ L 2d/(d-l) 

where we also have used that || • ||x,j> < || • \\ LP ' for p < p'. Now the threshold 
assumption (|3"2"|) and (fT?) show the boundedness of || w l h \\L°° and \\Vw l h \\ L 2d/(d-i). 

Finally, the proof of the case j = 1 in (|35l) is, with obvious changes, that of 
the equivalent result in J3 Lemma 3.1]. □ 

In the sequel we consider the auxiliary function vt ■ [0,T] — > V/, solution of 

v h + A h v h + U h F(u, u) = U h f, v h (0) = U h u . (37) 

According to jTCl Remark 4.2] we have 

uuw \\v h (t) - U h u(t)\\ < C\ \og{h)\h 2 , (38) 

for some constant C = C(M2). The following lemma provides a superconver- 
gence result. 

Lemma 2 Let (u,p) be the solution of (Qp-|0). Then, there exists a positive 
constant C such that the solution Vh of J 37\ l and the Galerkin approximation Uh 
satisfy the following bound, 

IMt)-«fc(t)||i < C\\og(h)\ 2 h\ te(0,T]. (39) 

1 /2 

Proof Since for = {vh — Uh) we have 

Vh + A h y h + A)l 2 Ii h (F{v h ,v h ) - F(u h ,u h )) = A 1 ^ 2 p h , 
where ph = Hh(F(vh, Vh) — F(u, u)), it follows that 

iMi)ii < rii4 /2e ~ (t ^iioii n ^ F ^'^)-^K,^))iio 

Jo 

+ f\\A h e-^ A -{A- h 1/2 Ph { s ))\\ ds. 
Jo 
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Applying ([33]) we have \\U h (F(v h ,v h ) ~ F(u h ,u h ))\\ < C \\yh\\ , so that taking 
into account that 

|| A l/2 e - (t - s )A h || o < (2e( ^ _ s ))-l/2 j ( 40 ) 

it follows that 

IllfcWIlo < 4= f^pk+ [)\A h e-^(^ 1/2 p h (s))\\ ds. 
\J le Jo \Jt — s Jo 

Since applying |161 Lemma 4.2] we obtain 

\A h e-^ A HA- 1/2 p h (s))\\ d S < C\\og(h)\ max \\ Ph ( S )\\ , 

u 0<s<t 



f 

Jo 



a generalized Gronwall lemma (26j pp. 188-189], together with (|33j) allow us to 
conclude 

IK - u h \\i < C\ \og(h)\\\v h - u\\o- 

Then by writing \\vh - u||o < \\vh - n^u||o + \\~HhU - u\\ and applying ([HJ 
and (f38"| , the proof is finished if we check that the threshold condition ([321 holds 
for w\ = Uh and w\ = Vh- In view of ([38| . (fT4[) and the inverse inequality (fl"2j) 
we have indeed that — u|| ■ = o(/i 3//2_J ), for j = 0, 1. In the case of u/j the 
threshold condition holds due to (l24l). □ 



Lemma 3 Let (w,j?) be the solution of ((I])-©. Then there exists a positive 
constant C such that 

\\v h (t)-u h (t)\\^<C\log{h)\ 2 h 2 , t€(0,T], (41) 

where Vh and Uh are defined by (|37|) and (|22p ~ (123|) respectively. 

Proof The difference f/j—M/j satisfies that Vh — Uh = Ah(vh— u/ l )+II^(£?(u, it) — 
B(uh,Uh), so that multiplying by A ft and taking norms, thanks to (|34l) . we 
have 

\\A^ 1/2 (v h - u h )\\ < \\Al /2 (v h - u h )\\ + C\\u - u h \\o ■ 
Now we write 

II" - u h \\ < \\u - U h u\\ Q + \\U h u - v h \\ Q + \\v h - u h \\ , 

so that (HU), (J3H1) and (|39j) , allow us to write, 

\\A- 1/2 (v h -u h )\\ <C\\og(h)\ 2 h 2 . 

Then, applying (f3Tj) the proof is finished. □ 

Lemma 4 Let (u,p) be the solution of ([l}-©. Then there exists a positive 
constant C such that 

\\^-i h m-i<^jjih r \]og(h)\ r ' , te(0,n r = 2,3,4, (42) 

where r' = 2 w/ien r = 2 and r' = 1 otherwise. 
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Proof The case r = 3,4 is proved in [THl Lemma 5.1]. For the case r = 2 we 
write 

u t - Uh = (u t - n h ti t ) + (U h u t - v h ) + (ih - u h )- (43) 

A simple duality argument and the fact that \\ut — HhUt\\o < C7j|| itt ||i, easily 
show that 

- nMi^ < ch^iutW, < c^h 2 . 

The bound of the third term on the right-hand side of (|43|) is given in Lemma [3l 
so that, thanks to the equivalence pip we are left with estimating 

y h =t 1 / 2 A- 1/2 (U h u t -v h ). 

We notice that 

Vh + A h y h = t l / 2 A)l 2 e h + ^r^ 2 A- 1/2 (n h u t - v h ), 
where Oh — (Ibj — Sh)u. Thus, 

Vh® = f s- x / 2 Al /2 e-^ A "(se h )ds 
Jo 

+ 1 f s- X/2 A)! 2 e-^ A »A- h \n h u s - v h ) ds. 
z Jo 

Recalling (|40|) by means of the change of variables r = s/t it is easy to show 
that 

£ s -./* K /V<-Mi^<^B(I,i), (44, 
where B is the Beta function (see e. g., [H]). Thus, we have 

\\Vhh < C B (\> l\™f* ( S INI° + \\A^(n h u s - v h )\\ ) . 

The first term on the right-hand side above is bounded by CM^h 2 . For the 
second one we notice hat 

A^ 1 {U h u t - v h ) = Oh - (n^u - v h ) 

so that using ([II ]) , (fT7 |) and i|ggj) it is bounded by M 2 h 2 |log(/i)|. □ 

Theorem 1 Let (u,p) be the solution of (QP-0). Then, there exists a positive 
constant C such that the postprocessed velocity u, defined in (126p - (|27[) . satisfies 
the following bounds: 

(i) Ifr = 2 then 

IKO-fi(**)lli<^^^|logW| 2 - (45) 

(ii) If r = 3,4 then 

IKO - u(t*)h < -7^^ r+1 - j |logWI, 3 = 0, 1. (46) 
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Proof The proof follows the same steps as [TBI Theorem 5.2]. Subtracting (|26[) 
from (TTJ, standard duality arguments show that 

< C(\\F(u(t*),u(t*))-F(u h (t*),u h (t*j)\\^ 1 + \\u t (t*)-u h (t*)\\^ 1 ). 

To bound the second term on the right-hand side above we apply Lemma 31 
whereas for the second we apply ([33]) to get 

\\F( U (t*), U (t*)) - F(u h (t*),u h (t*))\\-i < C\\u(t*) - u h (t*)h- (47) 

so that applying ([24)) the proof of (|45l) and the case j = 1 of (|46|) are finished. 
We now get the error bounds in the L 2 norm. It is easy to see that 

A(u(t*) - u(t*)) = U(F(u(t*),u(t*)) - F(u h (t*),u h (t*))) + U(ut(t*) - u h (t*)). 

Then, by applying A^ 1 to both sides of the above equations, we obtain 

||fi(f*) - u(t*)||o ^WA-iniFiuit*),^?)) - F(u h (t*),u h (t*)))\\ 
+ WA^Uiutit*) - u h (t*))\\ . 

As regards the nonlinear term, applying j fij Lemma 4.1] we obtain 

\\A^ 1 U(F(u(t*), u(f)) - F(u h (t*), u h (t*)))\\ < 

C(\\u(t*) - u h (i*)||_i + ||u(t*) - u h (t*)||i|K**) - u h (t*)\\ ). 

To bound the second term on the right-hand side above we apply (|24| , whereas 
the first one is bounded in the proof of [T|)J Theorem 5.2] by 

\\u(t*)-u h (t*)\\-i < ^)7i^ +1 |iog(MI- 

Finally, to bound ||^4 _1 II(u t (i*) — Uh(t*))\\o we apply [TBI Lemma 5.1] to obtain 

p-inMn - u h (r))||o < ^^i)7i^ +1 l iog(ft)|, 

which concludes the proof. □ 
In the following theorem we obtain the error bounds for the pressure p. 

Theorem 2 Let (u,p) be the solution of (Qp-flp. Then, there exists a positive 
constant C such that the postprocessed pressure, p, satisfies the following bounds: 

\\p(t*) -p(t*)hym < ^7/2^1 logWl r ' , (48) 
where r' = 2 if r = 2 and r' — 1 if r = 3, 4. 
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Proof The proof follows the same steps as [TSJ Theorem 5.3]. Applying the 
inf-sup condition ([5]) it is easy to see that 

f3\\p(t*)-p(t*)\\ L y M < ||«(f)-«(t*)||i + ||ut(<*)-A/.(«*)ll-i 

+ \\F(u h (t*),u h (t*)) ~ F(u(t*),u(t*))\\-i- 

Applying now (|45[) and (|46p to bound the first term and reasoning as in the 
proof of Theorem [T] to bound the other two terms we conclude (|48l) . □ 

Remark 1 As a consequence of Theorems [1] and [2] we obtain in the proof 
of Theorem [3] that (u — Uh) is an asymptotically exact estimator of the error 
(it — Uh) while (p — ph) is an asymptotically exact estimator of the error (p—ph)- 
However, as we have already observed u and p are not computable in practice. 
In Theorems [31 4 and [5] we present different procedures to get computable error 
estimators. 

As we pointed out before the MFE approximations (uh,Ph) to the velocity 
and the pressure of the solution (u,p) of the evolutionary Navier-Stokes equa- 
tions CO)-© a ^ an y nxe d time t* are also the approximations to the velocity 
and pressure of the steady Stokes problem P5)) - (P71) . In Theorem [3] we show 
that any a posteriori error estimator of the error in the steady Stokes problem 
(|26[) - (|27p gives us an a posteriori indicator of the error in the approximations 
to the evolutionary Navier-Stokes equations. 

Using the notation of [33] we will denote in the sequel by £ ve i({uh,Ph)j f, 
j = 0, 1, any a posteriori error estimator of the error Uh—u in the norm of IP(f2) d 
in the approximation to the velocity in the steady Stokes problem (|26p - (|27j) . 
We will denote by £, pre s{(uh,Ph), f, L 2 /M.) any error estimator of the quantity 

\\ph -pIUvr- 

Theorem 3 Let (u,p) be the solution o/([I])-([2]) and fix any positive timet* > 0. 
Assume that the Galerkin approximation (uh,Ph) satisfies, for h small enough 
and r — 2, 3, 4, 

|K**) - Uhit*)^ > C r h r -i, j =0,1. (49) 

for some positive constant C r — C r (t*). 

ft) V £vei((uh{t*),Ph{t*)), f-, H' J ), j — 0,1, is an efficient error indicator of 
the error in the MFE approximation to the steady Stokes problem (|26p - (127l) . 
That is, if there exist positive constants, C\ and C%, that are independent of the 
mesh size h, such that the following bound holds 

\\u(t*) - Uhit*)^ — V ' 

then t; ve i((uh(t*),ph{t*)), /, H 3 ), j = 0,1, it is also an efficient error indicator 
of the error in the MFE approximation to the evolutionary Navier-Stokes equa- 
tions, i.e. there exist positive constants C3 and C4 that are independent of the 
mesh size h such that the following bound holds for h small enough 

63 - ii — 77T\ 777771 < J=0,1. (51) 

|| u(t*) -U h (t*)\\j 
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(ii) If ^ ve i((uft(t*),pft(t*)), /, H 3 ), j = 0, 1 is an asymptotically exact error 
estimator of the error in the steady Stokes problem then it is also an asymp- 
totically exact error estimator of the error in the evolutionary Navier-Stokes 
equations. 

(Hi) Analogous results are obtained in the approximations to the pressure. 
In the case r — 2, the results are valid only in the H 1 norm. 

Proof For simplicity in the exposition we will concentrate on the cases r = 3, 4 
in the approximations to the velocity, the proof for the approximations to the 
pressure and for the case r = 2 being the same except for obvious changes. 
Let us first observe that 

IMP) - u(t*)\\ s < \\u h (t*) - u(t*)||j + IKO - u{t*)\\j, j = o, i. 

Dividing by \\uh(t*) — u(t*)\\j, using (j49|) and applying Theorem [T] we obtain 

, \\u h (t*)~u(t*)\\.j C t*- {{r - 1)/2 \^ 

Now, using (f50| we get 

\\u h (t*) - u(t*)\\j = IMP) - fi(g)Jjj gvoi(MP),Ph(P)), /, Hi) 
K(P) - IMP) - tv«\{{uh{t*),p h {t*))J,W) 

1 Z ve i((u h (t*), Ph (t*))J,W) 



< 



Ci \\u h (t*) - u(t*)\\j 



Taking h small enough so that ce l £ - h\log(h)\ < 1/2, we get 

2 - K(t*)- ' 1 j 

Now, we use the decomposition 

- "(P)||j < K(P) - u(i*)lli + IK**) - J = 0, 1. (53) 

Reasoning as before we get 

IM*->-*(OII, < 1 + ^^,1 log^i. 



K(t*) - «(t*)iij - a 

Since 

K(**)-«(f)llj > 1 Cvel(K(P),p h (P)),/,^') 



-«(f JHi - cr 2 ||«h(f)-«(i*)|| 3 - 

we finally reach 

£vei(K(OMP)),/,g J ) < 3^ ( , 4) 

From (JSll) and O we conclude flST]) with C 3 = Ci/2 and C 4 = 3C 2 /2. 
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Let us now assume that £ ve \((uh(t*),Ph(t*)), f, W) is an asymptotically ex- 
act error estimator. Using again the decomposition (|53[) we have 



K(o-fi(oii J =1 + lim pPf^pk = 1, 

the last equality being a consequence of Theorem (JTJ and the saturation hy- 
pothesis (|49p. As we pointed out before, this limit implies that (u — Uh) is an 
asymptotically exact estimator of the error (u — u^). Then 



h^o \\u h (t*)-u(t*)\\, 



= £ vc i((u h (t*), Ph {t*))J,Hi) |MQ -fi(**)llj = 

and t; ve \{{uh{t*) ,Ph{t*)) , f , Hi) is also an asymptotically exact estimator of the 
error in the approximation to the velocity of the evolutionary Navier- Stokes 
equations. □ 

Remark 2 We remark that with hypothesis (|49p we are merely assuming that 
the term of order h r ~i is really present in the asymptotic expansion of the 
Galerkin error. Let us also notice that the constant C r in (|49|) is, in general 
0(t*~ (r ~ 2)/2 ), so that the ratio £*~ ((r_1)/2) /Cr in the proof of Theorem 1 is, in 
general, 0{t* ( ' 1/2) ). 

The key point in Theorem[3]comes from the observation that if we decompose 

u — Uh = (u — u) + (u — Uh), (55) 

the first term on the right hand side of (|55[) . u — u, is in general smaller, by a 
factor of size 0(h \og(h)), than the second one, u — u% (Theorem [T]). Then, to 
estimate the error u — Uh we can safely omit the term u — u in (1551) . Compar- 
ing with the analysis of [33 for a nonstationary linear Stokes model problem 
the main difference is that the two terms in f|55[) are taken into account. In 
Theorem [4] below we show that this kind of technique can also be applied to 
the nonlinear Navier-Stokes equations. The advantage of this point of view is 
that hypothesis (|49"1) is not required for the proof of Theorem |U Let us finally 
observe that {uh,Ph) are the MFE approximations to the solution (u t ,Pt) of the 
Stokes problem that we obtain deriving respect to the time variable the Stokes 
problem (|26 ]> - ([27 |) ■ Then, we will denote by £ ve i((u h ,p h ), f t , H J), j = -1,0,1, 
any a posteriori error estimator of the error — u t in the norm of -£P (Q) d in the 
approximation to the velocity of the corresponding steady Stokes problem. The 
proof of the following theorem follows the steps of the proof of [5D1 Theorem 1] . 



Theorem 4 Let (u,p) be the solution of and let (uh,Ph) be Us MFE 

Galerkin approximation. Then, the following a posteriori error bound holds for 
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< t < T and a constant C independent of h. 

II (u - u h )(t)\\ < C\\u - u h (0)\\o + C£ vel ((u h (0),p h (0)), /(O), L 2 ) 
H vc i{{Mt),Ph{t)),.f(t),L 2 ) + Ct^ 2 max £ vel ((u h ,p h ), f, L 2 ) (56) 

0<s<t 

+Ct^ 2 max ^((ufc.Pfc),/.,^" 1 ). 

0<s<t 

Proof Let us denote by r\ — u — u. From ([2B)) - (f!?7]) it follows that 

rjt + Arj + IL(F(u, u) - F(u h , u h )) = U(u h - u t ). 
Then r\ satisfies the equation 

r,(t) = e- At r]{0) + f e~ A{t - s) IL(F{u, u) - F(u, u)) ds 
Jo 

+ [ e~ A ^~ s ^U{F{u hl u h ) - F{u,u) ds)+ [ e^^II^ - u t ) ds. 
Jo Jo 

Taking into account (J7]) we get 

„„<„„„ < ||, (0) || 0+c /' i^asp^* dt 

Jo Vt - s 

+c r \\A-^ 2 U(F(u h ,u h ) - F(u,u))\\o ds + c f p-^nK-uQUo ^ 

Jo y/t-S Jo y/t-S 

We first observe that for any v € L 2 {Vt) d we have || 17,2 ITf7 1 1 < C||«||_i. Then, 
taking into account (pJ5"j) we get 

\\A-V 2 IL(F{u,u) - F(u,u))\\ < C||fi-u|| 0l 
\\A-^ 2 U(F(u h ,u h )^F(u,u))\\ < C||u h -u|| . 

Let us observe that in order to apply ([55)1 we require uu to satisfy which 
holds due to (f24|). and ||m||oo and || Vu\\ L 2d/(d-i) to be bounded. Using ((4]) both 
norms are bounded in terms of \\u\\ 2. Applying (|19[) we get 

||w|| 2 < C(||uh,||o + \\u h ■ Vuft||o) 

< C (\\A h u h \\ + \\U h F(u h , u h )\\ + \\U h f\\ + \\u h ■ Vu h \\ ) . 

Finally, using that ||A^u/j||o is uniformly bounded, see (1281) . and reasoning as in 
(|29[) to bound the second and forth terms above we conclude \\u\\2 is uniformly 
bounded. Then, we arrive at 

iiiwiio < iw>)ii, +c f^kds+cf ds 

Jo \/t - s Jo \ft - s 

+C f ds. 
y/t~s 
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And then 



h(t)h < \\v(0)h + C [ ll^kds + Ct 1 / 2 max$ vel ((u h)Ph )J,L 2 ) 

Jo yt — s o<s<t 

+Ct 1/2 max €vei((u h ,Ph)J s ,L 2 ), 

0<s<t 

A standard application of a generalized Gronwall lemma [26 gives 
\\v(t)\\o < C|k(0)||o + C* 1/2 max £ ve i((u h ,p h )J, L 2 ) 

0<s<t 

+Ct 1/2 max £ ve i((uh,Ph),fs,L 2 ). 

0<s<t 

Now, using decomposition (|55[) we conclude the proof. □ 

We observe that using the same proof, a similar bound for the norm of 

the error can be obtained changing only £ ve] ( (u h ,p h ), f, L 2 ) by £ vcl ( (u h , p h ) , /, H 1 ) 
and £ ve \((uh,Ph),ft, H^)hy £ vc i((uh,Ph), ft, L 2 ). Let us also remark that The- 
orem 2] allows to a posteriori obtain upper error bounds for the error in the 
approximation to the nonlinear Navier-Stokes equations using only upper error 
bounds for some Stokes problems depending only on the data and the com- 
puted approximation. However, the estimation of the error at a time t requires 
the estimation of the error of a family of Stokes problems with right hand side 
depending on r, for all r G [0, i]. 

We now propose a simple procedure to estimate the error which is based on 
computing a MFE approximation to the solution (u(t*),p(t*)) of (|26|) - ([27f on a 
MFE space with better approximation capabilities than (X^ jT ., Qh,r-i) in which 
the Galerkin approximation (uh,Ph) is defined. This procedure was applied to 
the p- version of the finite-element method for evolutionary convection-reaction- 
diffusion equations in |19j . The main idea here is to use a second approximation 
of different accuracy than that of the Galerkin approximation of (u, p) and whose 
computational cost hardly adds to that of the Galerkin approximation itself. 

Let us fix any time t* G (0, T] and let approximate the solution (u,p) of the 
Stokes problem (|26]) - ((27l) by solving the following discrete Stokes problem: find 
(u h (t*),ph(t*)) G X x Q satisfying 

(v« h (t*),V0) + (vji h (O,0) - (f,i)-(F(u h (t*),u h (t*)),fy (57) 

(v-« fc (f),^) = V-0 G Q, (58) 

where (X, Q) is either: 
(a) The same-order MFE over a finer grid. That is, for h! < h, we choose 

(X,Q) — {Xh\riQh',r-l)- 
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(b) A higher-order MFE over the same grid. In this case we choose (X, 

Q) = {Xh,r+l,Qh,r)- 

We now study the errors u — Uh and p — Ph- 

Theorem 5 Let (u,p) be the solution of ([I])-© an d f or r — 2,3,4, and let 
(|10p - (|lip hold with k = r + 2 Then, there exists a positive constant C such that 
the postprocessed MFE approximation to u, Uh satisfies the following bounds for 
r = 2,3,4 and t € (0,T]: 

(i) if the postprocessing element is (X,Q) = (Xh', r ,Qh>,r-i), then 

\\u(t) - Mm, < ^Sm {h ' r ^ + ijrm hr+1 ~ 3 \ lo s W r '' j = °> ^ 

(59) 

im-PhmLyu^^^ihr-' + ^^h^iogihr', (eo) 

(ii) if the postprocessing element is (X,Q) — (Xh, r +i,Qh,r), then 

\\u(t) - u^t)^ < J( ^ 2 h^\\og{h)f 1 j =0,1, (61) 

\\p(t)-p h (t)\\ L 2 /R < -^t^H log(h)| r '. (62) 

For r — 2 only the case j — 1 in i59\) and i61\) holds. In l59\) - ![b!^) . r' = 2 when 
r = 2 and r' = 1 otherwise. 

Proof The cases r = 3,4 have been proven in Theorems 5.2 and 5.3 in [T6] . 
Following the same arguments, we now prove the results corresponding to r = 2 
and (X, Q) = (Xh>,r>Qh',r-i)> the case (X, Q) — (X h ^ r+ i,Q htr ) being similar, 
yet easier. We decompose the error u — Uh = (u — s^) + (s^ — %), where 
(sh')Qh') € x (3,,/ :1 is the solution of 

(Vs h /,V<^)-(<Zh',V-<^) = (/-^.uj-ut^v) V0v GJCv,a,(63) 
(V-s/^Vv) - V^eQw, (64) 

that is, Sh' is the Stokes projection of u onto V^/. Since in view of ([T7 ]) -(fT8 |) we 
have 

ll u - s h'\\i + \\p- qh>\\ L 2 /R < CM 2 ti, 

we only have to estimate sw — Uh and qw — Ph- To do this, we subtract (|5T[) 
from l|63p . and take inner product with e-h = sv — Uh to get 

||Ve/i|| < (\\u t - ii h \\_ 1 + \\F(u h ,u h ) - F(k,u)||_ 1 ) ||e/j.|| x . 

Now applying Lemma [U f|35[) and (f24| the proof of (f59|) is finished. 

To prove ([60]) . again we subtract (|57|) from |63|) . rearrange terms and apply 
the inf-sup condition ([5]) to get 

P\\lh> -Ph\\L 2 /R < HVe/iHo + ||uf - UhlUx + \\F(u h ,u h ) - F(u,u)||_ 1 

and the proof is finished with the same arguments used to prove (|59|) . □ 
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To estimate the error in (uh{t*),ph(t*)) we propose to take the difference 
between the postprocessed and the Galerkin approximations: 

fjh,vei(t*) = u h (t*) - u h (t*), fjh, prea (t*) = Ph(t*) -Ph(t*)- 

In the following theorem we prove that this error estimator is efficient and 
asymptotically exact both in the L 2 {Vl) d and norms and it has the 

advantage of providing an improved approximation when added to the Galerkin 
MFE approximation. 

Theorem 6 Let(u,p) be the solution of [Ity-^ and fix any positive timet* > 0. 
Assume that condition J4ff| ] * s satisfied. Then, there exist positive constants ho, 
70 < 1, and C\, C2, C3 and C4 such that, for h < hg and < 7 < 70, the 
error estimators fjh,vei(t*) fjh,pies(t*) satisfy the following bounds when (X, Q) = 
(Xh', r ,Qh',r-i) and hi < 7/1: 



Ci < 



\\Vh,vc\(t*)\\j 

l(tt-tt/i)(**)llj 



<C 2 , j = 0,1, c 3 < 



|%,pres(OI!-L 2 /R 



< C 4 . (65) 



Furthermore, if(X,Q) = {Xh', r ,Qh',r-i) 
(X h ,r+i,Qh,r) then 



\\(p-Ph)(t*)\\ L 2 
with h! = h 1+e , e > 0, or (X,Q) 



lim 



|%,vel(i*)||. 



= 1, j = 0, 1, lim 



\fjh,pi-Qs(t*)\\L 2 /R 



= 1. 



h^o \\{u - u h ){t*)\\j ,J '' fc-+o \\(p-p h )(t*)hy» 
For the mini element, the case j = in 165\) and i66\) must be excluded. 



(66) 



Proof We will prove the estimates for the velocity in the case r = 3, 4, since the 
estimates for the pressure and the case r — 2 are obtained by similar arguments 
but with obvious changes. Let us observe that for j = 0, 1 



HO -MOIL, 



< 
< 



|%,vel(t*)||j + ||«h(**)-w(**)||j 



+ 



c 



(t*)(r-i)/: 



(£*)(r-2)/2 

7 h r+1 - j \\og(h)\. 



On the other hand 

||%,vel(** 



< 
< 



IKO-«h(**)lli + ll«/.(**)-«(Olli 

\\ U (t*)-u h (tn\\ j + 7 —^ m (hr- j 



+ 



c 



( t *)(r-l)/: 



(f*)(r-2)/2 

r /T +1 ^|log(/i)|. 



Using (149)) we get 

Ifevel (**)||,' 



|(tl - U h ){t* 



< 



-(r-2)/2 I ^ T 



(t*)- (r - 1)/2 \log(h)\h 



(67) 
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Taking h' < 7/1 and h and 7 sufficiently small, the bound is readily obtained. 
The proof of (|66|) follows straightforwardly from (|67|) . since in the case when 
(X, Q) = (X h ', r , Qh',r-i) with ti = h 1+e , e > 0, the term (ti /h) r ~ j when h 
tends to zero, and in the case when (X, Q) = (X^ , r +i, Qh.r) the term containing 
the parameter ti is not present. □ 



4 A posteriori error estimations. Fully discrete 
case 

In practice, it is not possible to compute the MFE approximation exactly, and, 
instead, some time-stepping procedure must be used to approximate the solu- 
tion of ([22)- (EH- Hence, for some time levels = t < h < ■ ■ ■ < t N = T, 
approximations UJf ~ Uh(t n ) and P£ ~ Ph(t n ) are obtained. In this section we 
assume that the approximations are obtained with the backward Euler method 
or the two-step BDF which we now describe. For simplicity, we consider only 
constant stepsizes, that is, for N > 2 integer, we fix k — T/N, and we denote 
t n = nk, n = 0, 1, . . . , N. For a sequence (y n )n=o we denote 

Dy n = y n -y n -\ n=l,2...,N. 

Given = u/i(0), a sequence (U£,P£) of approximations to (uh(t n ),Ph(t n )), 
n = 1, . . . N, is obtained by means of the following recurrence relation: 

(d t [/£,^) + (VC/,?,V^) (68) 

+ b (uz, u%, 4> h ) - (if, v ■ M = (/, <M v<f> h e x h . r , 

(V-E#,^) = 0, V^ h e Qfc ir -i, (69) 

where d t — k~ 1 D in the case of the backward Euler method and d t = k~ 1 (D + 
\D 2 ) for the two-step BDF. In this last second starting value U\ is 

needed. Here, we will always assume that U\ is obtained by one step of the 
backward Euler method. Also, for both the backward Euler and the two-step 
BDF, we assume that U® = Uh(0), which is usually the case in practical situa- 
tions. 

We now define the time-discrete postprocessed approximation. Given an 
approximation d^U]^ to Uh(t n ), the time-discrete postprocessed velocity and 
pressure (U n , P n ) are defined as the solution of the following Stokes problem: 

(yu n , V4>) + (vp'\ 0) = (/, 0) - b (u%, un,4>) - (d*uj;,ct>) , e ff W, 

(70) 

(v • u n ,ip) = 0, w e L 2 (n)/R. (71) 

For reasons already analyzed in [TS] and [T7] we define 

d\Ul = IL h f - A h Ul - K h F Ut) (72) 
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as an adequate approximation to the time derivative Uh (t n ) . 

For the analysis of the errors u{t) — U n and p(t) — P n we follow [I7j . where 
the MFE approximations to the Stokes problem (|70l) - ([TTj) are analyzed. We 
start by decomposing the errors u(t) — U n and p(t) — P" as follows, 

u(t n )-U n = (u(t)-u(t n )) + e n , (73) 
p(t n )-P n = (p(t n )-p(t n ))+n n , (74) 

where e" = {i(< n ) — t/™ and 7r™ = p{t n ) — P n arc the temporal errors of the time- 
discrete postprocessed velocity and pressure (U n ,P n ). The first terms on the 
right-hand sides of (I73l - ([74]l are the errors of the postprocessed approximation 
that were studied in the previous section. 

Let us denote by eJJ = Uh(t n ) — Uj£, the temporal error of the MFE approx- 
imation to the velocity, and by — ph{t n ) — the temporal error of the 
MFE approximation to the pressure. In the present section we bound (e™ — e^) 
and (7f n — 7r£) in terms of e£. 

The error bounds in the following lemma are similar to those of |17l Propo- 
sition 3.1] where error estimates for MFE approximations of the Stokes prob- 
lem CnH-CEE]) are obtained. 

Lemma 5 There exists a positive constant C = C(maxo<t<T ||^4/iUft(i)||o) such 
that 

W^-ell < Ch 2 -1(\\el\\ 1 + \\el\\l + \\A h el\\ ), j = 0,l, 1 < n <(72S) 

W^-n^hyR < CMKIIi + KII? + ||AeJ*||o), 1 < n < N. (76) 

Proof Let us denote by I = g + (d*[/™ - u h (t n )) where g = F(U%,U%) - 
F(u h (t n ),u h (t n j). Subtracting dZOJ) — (LTD) from we have that the tem- 

poral errors (e™, jt n ) of the time-discrete postprocessed velocity and pressure are 
the solution of the following Stokes problem 

(Ve",V</>) + (W",0) = (?,</>), V0€J3oW> (77) 
(V-e",V) = 0, V^I 2 (fl)/l. (78) 

On the other hand, subtracting (f68|) -(l6"9 | from (|22]) - ([23|) and taking into ac- 
count that, thanks to definition (|72l) dtU£ — d^UJ^, we get that the temporal 
errors (e^,7r^) of the fully discrete MFE approximation satisfy 

(Ve£, V4> h ) + (V7r£,<M) = (l,4> h ), V0 h G X hlT , 
(V-e^,Vh) = 0, Vip h e Q h , r -i, 

and thus (eJJ,7r^) is the MFE approximation to the solution (e n ,n n ) of ([77|- 
(|75|) . Using then dHJ) we get 

\\e--e^\ j <Ch 2 ^\\l\\ . 

For the pressure we apply (|18p and (f2"U| to obtain 

\\n n - Khym < Ch\\w n \\ m/R < Ch\\l\\ . 
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Then, to conclude, it only remains to bound ||i||o- From the definition of d* t XJ^ 
it is easy to see that 

d* t U% - u h (t n ) = A h el - U h (F ([/«, Eft) - F(u h (t n ), u h (t n ))) , 

so that 

|Kt/£-^(i n )|| < ll^Ho + lbllo. 

Now, by writing g as 

g = F(e n h , u h (t n )) + F(u h (t n ), e n h ) - F(e n h , eft, 
and using ([2"9" ]) -([3"0 ]) we get 

llffllo < (P^(t„)||o||4||i + ||4||? /2 p, ie ^||y 2 ) , 

from which we finally conclude (|75p and ([To]) . □ 

Let us consider the quantities U n — U£ and P" — P£ as a posteriori indicators 
of the error in the fully discrete approximations to the velocity and pressure 
respectively. Then, we obtain the following result: 

Theorem 7 Let (u,p) be the solution of and let (0J) hold. Assume 

that the fully discrete MFE approximations {U^, P^), n = . . . , N = T/k are 
obtained by the backward Euler method or the two-step BDF i68\) -![b^ \) . and let 
(U n , P") be the solution of (7§-(7^l. Then, for n = 1, . . . , N, 

\\U n -Ul% < \\u(t n )-u h (t n )\\ 3 +C' la h 2 '^, j — 0, 1, (79) 

\\P n -Pn L yM < \\p{tn)-p h {tn)\\LVR + C[ h^, (80) 

tn 

where C[ o is the constant in = 1 for the backward Euler method 

and Iq — 2 for the two-step BDF. 

Proof In 17, Theorems 5.4 and 5.7] we prove that if (|9]) and the case I = 2 
in ([17]) hold, the errors eJJ of these two time integration procedures satisfy for 
k small enough that 

||e£||o + U^I|o<C*o7^ l<n<AT, (81) 

for a certain constants C\ and C2, where Iq = 1 for the backward Euler method 
and l = 2 for the two-step BDF. Since ||A^ /2 e£|| < ||eJ>||S /2 ||A,,e]>||£ /2 , and 
then ||e£||i < C\\el\\l /2 \\A h el\\l /2 , from flU]) and (ffB) - (^6|l we finally reach that 
for k small enough 

\\~e n - el\\, < Ciy-'^, j = 0, 1, 1 < n < N, (82) 
W^-Khyw < C h h^, l<n<N, (83) 
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where C[ o is a positive constant. 







Let us decompose the estimators as follows: 




(84) 



(85) 



yr-UZWi < \\u{t n )-u h {t n )\\ j + \\el-e n \\ h j = 1,2, 
||P" - PZhyR < ||p(*») ~ Ph(t n )\\ L y R + - n n \\ L */*- 



Thus, in view of ([52" |) -(|55 ]) we obtain JTS} and §U§ 



□ 



Let us comment on the practical implications of this theorem. Observe 
that from ([84]) and flU) the fully discrete estimators U n - U% and P™ - P£ 
can be both decomposed as the sum of two terms. The first one is the semi- 
discrete a posteriori error estimator we have studied in the previous section (see 
Remark [lj and which we showed it is an asymptotically exact estimator of the 
spatial error of U£ and P£ respectively. On the other hand, as shown in (|52"j)- 
(|83[) . the size of the second term is in asymptotically smaller than the temporal 
error of UJl and P% respectively. We conclude that, as long as the spatial an 
temporal errors are not too unbalanced (i.e., they are not of very different sizes), 
the first term in and (|55|) is dominant and then the quantities U n — U£ 
and P n — P£ are a posteriori error estimators of the spatial error of the fully 
discrete approximations to the velocity and pressure respectively. The control 
of the temporal error can be then accomplished by standard and well-stablished 
techniques in the field of numerical integration of ordinary differential equations. 

Now, we remark that (U n ,P n ) are obviously not computable. However, 
we observe that the fully discrete approximation {U£,P£) of the evolutionary 
Navier-Stokes equation is also the approximation to the Stokes problem (1701) - 
(fTTj) whose solution is (U n ,P n ). Then, one can use any of the available error 
estimators for a steady Stokes problem to estimate the quantities \\U n — U%\\j 
and || P" — P^\\l 2 /ffi> which, as we have already proved, are error indicators of the 
spatial errors of the fully discrete approximations to the velocity and pressure, 
respectively. 

To conclude, we show a procedure to get computable estimates of the error 
in the fully discrete approximations. We define the fully discrete postprocessed 



(w£, V0) + (vp£, Pj = (/, $) - b (un, uzJ) - (d* t uj;, 4>) e x, 





(86) 
(87) 
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where (X, Q) is as in (j5T|) - (l5g]) . 

Let us denote by e% = Uh{t n ) — U% and tt^ = Ph{t n ) — P\l the temporal 
errors of the fully discrete postprocessed approximation {U^, P£) (observe that 
the semi-discrete postprocessed approximation {uh,Ph) is defined in (|57|) - (|58l ). 
Let us denote, as before, by eJJ the temporal error of the MFE approximation 
to the velocity. Then, we have the following bounds. 

Lemma 6 There exists a positive constant C = C (m&xo<t<T \\AhUh(t)\\o) such 
that for 1 < n < N the following bounds hold 

\K-el\\ 3 < C^OKIIi + lle^lf + H^llo), j — 0, 1(88) 
\K-K\\ma)/M < Ch(\Wh + \M\\l + \\Ah<%\\o). (89) 

Proof The bound (|88]) is proved in [TTl Proposition 3.1]. To prove (189)) we 
decompose 

\\K - K IU 2 (n)/R < \\K - 7T n Hi 2 (n)/M + ||7r" - K\\L2(n)/R- 

The second term above is bounded in (|76|) of Lemma [5l For the first we observe 
that 7f^ is the MFE approximation in Q to the pressure Tr n in (f77]) - (f78l) so that 
the same reasoning used in the proof of Lemma [5] allow us to obtain 

R - ff n |U» ( n)/n < Ch\\n n \\ m/R < CfcflKUi + || e ,"J? + \\A h el\\ ). 

a 

Using (181) as before, we get the analogous to (|82|) and (|83|) . i.e., for k small 
enough the following bound holds 

\\e n h - eft\j < C' l0 h 2 -^, j = 0, 1, 1 < n < N, (90) 

WK-Khym < C' lo h^, l<n<N. (91) 

tn 

where C[ Q is a positive constant. 

Similarly to J83-(g5) we write E7£ - U£ = (u h (t n ) - u h {t n )) + (eg - gjj) 

and P™ - P™ = (ph(t n ) ~ Ph(t n )) + (K - TTfc), so that in view of ([gO ]) -(|9"T j) wc 
have the following result. 

Theorem 8 Let (u,p) be the solution of and let (0) hold. Assume 

that the fully discrete MFE approximations (U£, Pfi), n = . . . , N = T/k are 
obtained by the backward Euler method or the two-step BDF l6'ff)) -i l7>PI) , and let 
(U£,P%) be the solution of fSfy-gflj). Then, for n = 1, . . . , N, 

WVZ-UJtWi < \\uh(tn)-u h (t n )\\ j+ C' lo h 2 -^, i = 0,l, (92) 

\\Ph-Phhyn < \\Mtn)~p h (t n )\\Lym + C" lo h^, (93) 
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where C' lg is the constant in i90\) - [§!]) , Iq = 1 for the backward Euler method 
and Iq — 2 for the two-step BDF. 

The practical implications of this result are similar to those of Theorem 
that is, the first term on the right-hand side of (j9"2")l is an error indicator of the 
spatial error (see Theorem [6J while the second one is asymptotically smaller 
than the temporal error. As a consequence, the quantity (U% — U%) is a com- 
putable estimator of the spatial error of the fully discrete velocity U£ whenever 
the temporal and spatial errors of U£ are more or less of the same size. As be- 
fore, similar arguments apply for the pressure. We remark that having balanced 
spatial and temporal errors in the fully discrete approximation is the more com- 
mon case in practical computations since one usually looks for a final solution 
with small total error. 

As in the semi-discrete case, the advantage of these error estimators is that 
they produce enhanced (in space) approximations when they are added to the 
Galerkin MFE approximations. 

5 Numerical experiments 

We consider the equations 

u t - vAu + (u • V)li + Vp = /, (94) 
div(u) = 0, 

in the domain fl — [0, 1] x [0, 1] subject to homogeneous Dirichlet boundary 
conditions. For the numerical experiments of this section we approximate the 
equations using the mini-element jlOj over a regular triangulation of f2 induced 
by the set of nodes (i/N,j/N), < i,j < N, where N = 1/h is an integer. 
For the time integration we use the two-step BDF method with fixed time 
step. For the first step we apply the backward Euler method. In the first 
numerical experiment we study the semi-discrete in space case. To this end in 
the numerical experiments we integrate in time with a time-step small enough 
in order to have negligible temporal errors. We take the forcing term f(t,x) 
such that the solution of (IM1) with v = 0.05 is 

u 1 (x,y,t) = 2nip(t) sin 2 (7ra;) sin(7n/) cos(7ry), 
u 2 (x,y,t) = — 2Tnp(t) sin 2 (7ry) sin(7ra;) cos(7ra;), (95) 
p(x,y,t) = 2Q v (t)x 2 y. 

We chose ip(t) = t in the first numerical experiment. 

When using the mini-element it has been observed and reported in the lit- 
erature (see for instance [3D], IS]) H [31]) [SI] and [3S]) that the linear part of 
the approximation to the velocity, u l h , is a better approximation to the solution 
u than Uh itself. The bubble part of the approximation is only introduced for 
stability reasons and does not improve the approximation to the velocity and 
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pressure terms. For this reason in the numerical experiments of this section 
we only consider the errors in the linear approximation to the velocity. Also, 
following [5], we postprocess only the linear approximation to the velocity, i.e., 
we solve the Stokes problem (|S"7|) - (j5"51l with u l h and ii l h on the right-hand-side 
instead of Uh and iih- The finite element space at the postprocessed step is the 
same mini-element defined over a refined mesh of size h! . We show the Galerkin 
errors and the a posteriori error estimates obtained at time t* = 0.5 by taking 
the difference between the postprocessed and the standard approximations to 
the velocity and the pressure. In Figure [TJ we have represented the errors in the 




Figure 1: Errors (solid lines) and estimations (dashed lines) in L 2 (asterisks) 
and H 1 (circles) for h = 1/10, 1/12, 1/14, 1/16 and 1/18 and b! = 1/24, 1/30, 
1/34, 1/38 and 1/40 respectively. On the left, error estimations for the first 
component of the velocity. On the right, error estimations for the pressure. 

first component of the velocity of the Galerkin approximation in the L 2 and H 1 
norms and the errors for the pressure in the L 2 norm using solid lines. We have 
used dashed lines to represent the error estimations. The results for the second 
component of the velocity are completely analogous and they are not reported 
here. The L 2 errors of the pressure, on the right of Figure Q] are approximately 
twice as those of the H 1 errors of the velocity, on the left of Figure [IJ in this 
example. We can observe that with the procedure we propose in this paper 
we get very accurate estimations of the errors, specially in the H 1 norm of the 
velocity. The difference between the behavior of the error estimations in the 
L 2 and H 1 norms of the velocity are due to the fact that for first order ap- 
proximations the postprocessed procedure increases the rate of convergence of 
the standard method only in the H 1 norm for the velocity and the L 2 norm 
for the pressure. However, since the postprocessed method produces smaller 
errors than the Galerkin method also in the L 2 norm it can also be used to 
estimate the errors in this norm, as it can be checked in the experiment. On the 
right of Figure[T]we can clearly observe the asymptotically exact behavior of the 
estimator in the L 2 errors in the pressure in agreement with (|66|) of Theorem |6l 
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Let us denote by 



n _ ul(t*)-uUt*) _ Ph(t*) - p(t*) 

vel ~ ui(t*)-u{(t*)' ° pie ~ p(t*)- Ph (t*y 

the efficiency indexes for the first component of the velocity and for the pressure. 
In Table[T]we have represented the values of the L 2 and H 1 norms of the velocity 
index and the L 2 /R norm of the pressure index for the experiments in Figure [TJ 
We deduce again from the values of the efficiency indexes that the a posteriori 



h 


#vcl 


#vcl 1 


#prc L 2 /R 


1/10 


1.3640 


0.7721 


1.2588 


1/12 


1.3280 


1.0197 


1.1602 


1/14 


1.1695 


1.0068 


1.1084 


1/16 


1.3259 


0.9290 


1.0526 


1/18 


1.2741 


1.0438 


1.0167 



Table 1: Efficiency indexes 

error estimates are very accurate, all the values are remarkably close to 1, which 
is the optimal value for the efficiency index. More precisely, we can observe 
that the values of the efficiency index in the L 2 norm for the velocity in this 
experiment belong to the interval [1.1695, 1.3640]. The values in the H 1 norm 
for the velocity lie on the interval [0.7721, 1.0438] and, finally, the values for the 
pressure are in the interval [1.0167, 1.2588]. 

To conclude, we show a numerical experiment to check the behavior of the 
estimators in the fully discrete case. We choose the forcing term / such that 
the solution of ([94]) is ([95]) with ip(t) = sm((2n + n/2)t). The value of v = 0.05 
and the final time t* = 0.5 are the same as before. In Figure [2] on the left, we 
have represented the errors obtained using the implicit Euler method as a time 
integrator for different values of the fixed time step k ranging from k = 1/10 to 
k = 1/160 halving each time the value of k. For the spacial discretization we 
use the mini-element with always the same value of h = 1/18. We use solid lines 
for the errors in the Galerkin method and dashed lines for the estimations, as 
before. The L 2 norm errors are marked with asterisks while the H 1 norm errors 
are marked with circles. We estimate the errors using the postprocessed method 
computed with the same mini-element over a refined mesh of size h! — 1/40. 
We observe that the Galerkin errors decrease as k decreases until a value that 
corresponds to the spatial error of the approximation. On the contrary, the 
error estimations lie on an almost horizontal line, both for the velocity in the 
L 2 and H 1 norms and for the pressure. This means, as we stated in Section 
4.2, that the error estimations we propose are a measure of the spatial errors, 
even when the errors in the Galerkin method are polluted by errors coming 
from the temporal discretization. In this experiment the error estimations are 
very accurate for the spatial errors of the velocity in the H 1 norm and for the 
errors in the pressure. As commented above, the fact that postprocessing linear 
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Figure 2: Errors (solid lines) and estimations (dashed lines) in L 2 (asterisks) 
and H 1 (circles) for h = 1/18. On the left: Euler; on the right: two-step BDF 
for k = 1/10 to k = 1/160. 



elements does not increase the convergence rate in the L 2 norm is reflected in 
the precision of the error estimations in the L 2 norm. On the right of Figure [2] 
we have represented the errors obtained when we integrate in time with the 
two-step BDF and fixed time step. The only remarkable difference is that, as 
we expected from the second order rate of convergence of the method in time, 
the temporal errors are smaller for the same values of the fixed time step k. 
Again, the estimations lie on a horizontal line being essentially the same as in 
the experiment on the left. 
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